ASYMMETRIC DRIFT 

AND THE STELLAR VELOCITY ELLIPSOID 



Kyle B. Westfall\ Matthew A. Bershady^ Marc A. W. Verheijen^, 
David R. Andersen'^, & Rob A. Swaters"^ 

University of Wisconsin - Madison, Kapteyn Institute, 
'^NRC Herzherg Institute of Astrophysics, ^ University of Maryland 



Abstract We present the decomposition of the stellar velocity ellipsoid using stellar ve- 
locity dispersions within a 40° wedge about the major-axis (crmaj), the epicycle 
approximation, and the asymmetric drift equation. Thus, we employ no fitted 
forms for cTmaj and escape interpolation errors resulting from comparisons of the 
major and minor axes. We apply the theoretical construction of the method to 
integral field data taken for NGC 3949 and NGC 3982. We derive the vertical-to- 
radial velocity dispersion ratio (oz/an) and find (1) our decomposition method 
is accurate and reasonable, (2) NGC 3982 appears to be rather typical of an Sb 
type galaxy with a^/oR = O.TSIq j"^ despite its high surface brightness and 
small size, and (3) NGC 3949 has a hot disk with a, /an = l.lSto.lt- 



1. Motivation and Methodology 

The shape of the stellar velocity ellipsoid, defined by aji, a^, and Uz, pro- 
vides key insights into the dynamical state of a galactic disk: Uz'-crR provides a 
measure of disk heating and a^'.cTR yields a check on the validity of the epicy- 
cle approximation (EA). Additionally, ctr is a key component in measuring the 
stability criterion and in correcting rotation curves for asymmetric drift (AD), 
while (Tz is required for measuring the disk mass-to-light ratio. The latter is 
where the DiskMass survey focuses (Verheijen et al. 2004, 2005); however, in 
anything but face-on systems, cr^ must be extracted via decomposition of the 
line-of-sight (LOS) velocity dispersion. Below, we present such a decomposi- 
tion for two galaxies in the DiskMass sample: NGC 3949 and NGC 3982. 

Previous long-slit studies (e.g., Shapiro et al. 2003 and references therein) 
acquired observations along the major and minor axes and performed the de- 
composition via the EA and AD equations; using both dynamical equations 
overspecifies the problem such that AD is often used as a consistency check. 
Here, use of the SparsePak (Bershady et al. 2004, 2005) integral field unit 
(IFU) automatically provides multiple position angles, thereby increasing ob- 
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serving efficiency and ensuring signal extraction along the desired kinematic 
axes. Long-slit studies have also used functional forms to reduce the sensitiv- 
ity of the above decomposition method to noise. Here, only measures of the 
LOS velocity dispersions within a 40° wedge about the major axis are used 
to perform the decomposition by incorporating both the EA and AD equations 
under some simplifying assumptions. Velocities and radii within the wedge 
are projected onto the major axis according to derived disk inclinations, i, and 
assuming near circular motion. In the end, our method requires neither fitted 
forms nor error-prone interpolation between the major and minor axes. Future 
work will compare this decomposition method with the multi-axis long-sUt 
method and investigate effects due to use of points off the kinematic axes. 

Following derivations in Binney & Tremaine (1987) and assuming (1) EA 
holds, (2) the velocity ellipsoid shape and orientation is independent of z 
{d{vRV^)/dz = 0), (3) both the space density, u, and aji have an exponen- 
tial fall off radially with scale lengths of hn and 2hji, respectively, and (4) the 
circular velocity is well-represented by the gaseous velocity, Vg, the equation 

for the AD of the stars becomes Vg-lu'^ = ^aj^ + where 

is the mean stellar rotation velocity; hence, aji is the only unknown. The third 
assumption requires mass to follow light, S oc a^, and constant velocity ellip- 
soid axis ratios with radius; S is the surface density. The major-axis dispersion 
is geometrically given by cr^aj = a'^sin^i + {rjanj'^cos^i, where rj = az/o-R 

is constant with radius. Finally, EA, o"^/(t|. = | ( + l) ' completes a full 
set of equations for decomposition of the velocity elhpsoid. 

2. Analysis 

Data for testing of the above formalism was obtained during SparsePak 
commissioning (Bershady et al. 2005, see Table 1). Both NGC 3949 and 
NGC 3982 were observed for 3 x 2700s at one IFU position with = 513.1 
nm and A/AA = 11700 or 26 km s~^. The velocity distribution function 
(VDF) of both gas and stars is parameterized by a Gaussian function. In each 
fiber with sufficient signal-to-noise, the gaseous VDF is extracted using fits to 
the [Oiii] emission line; the stellar VDF is extracted using a modified cross- 
correlation method (Tonry & Davis 1979; Statler 1995) witii HR 7615 (KOIII) 
as the template (Westfall et al. 2005). The pointing of the IFU on the galaxy is 
determined post factum to better than 1" by minimizing the difference be- 
tween the fiber continuum flux and the surface brightness profile. Subsequent 
galactic coordinates have been deprojected according to the kinematic i and 
position angle. 

Figure I shows LOS dispersions for both the gas, ag^LOS, and stars, cr*,L05; 
data points are given across the full field of the IFU with points along the major 
and minor axes and in between having different symbols (see caption). From 
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this Figure note (1) there is no significant difference in glos along the major 

and minor axes for NGC 3949 and (2) the large gas dispersion within R < 10" 

for NGC 3982 is a result of poor single Gaussian fits to the multiple dynamical 

components of its LINER nucleus. 
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Figure 1. Measurements of the LOS dispersions in NGC 3949 (left) and NGC 3982 (right) 
for stars (top) and gas (bottom) as a function of galactic radius; the top axis is in units of scale 
length (hn is 20 and 10 arcsec in NGC 3949 and NGC 3982, respectively). Major axis, minor 
axis, and in between fibers are shown as points, triangles, and crosses, respectively; "axis fibers" 
lie in a 40° wedge around the axis. The binned mean values are shown for the major axis, minor 
axis, and otherwise as dashed, dotted, and solid lines, respectively. 





Figure 2 gives the folded gaseous and stellar rotation curves and compares 
the measured cTmaj from Figure 1 with that calculated using the formalism 
from §1. The value of rj used in Figure 2 provides the minimum difference 
between the two sets of data (as measured by x^; see Figure 3a). We find 
r? = and rj = O.JStoH for NGC 3949 and NGC 3982, respectively; 

errors are given by 68% confidence limits. A comparison of these values to 
the summary in Shapiro et al. (2003) is shown in Figure 3b. The disk of NGC 
3982 is similar to other Sb types studied; however, NGC 3949 seems to have 
an inordinately hot disk. The latter, while peculiar, is also supported by the 
indifference between its major and minor axis aLOS from Figure 1 and the 
same indifference seen in the Call data presented by Bershady et al. (2002). 
Our streamlined velocity-ellipsoid decomposition method appears accurate, as 
seen by comparison with (1) galaxies of a similar type for NGC 3982, and (2) 
data in a different spectral region for NGC 3949. 

This work was generously supported by an AAS International Travel Grant, 
the Wisconsin Space Grant Consortium, and the NSF Grant AST-0307417. 
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Figure 2. Comparison of measured and calculated amaj in NGC 3949 (left) and NGC 3982 
(right). Top panels show the deprojected rotational velocity for both the gas (grey triangles) and 
stars (Mack points) along the major axis. Bottom panels compare the measured amaj (points) 
to those from our formalism using rj as derived from Figure 3a (squares; offset in R for compar- 
ison). Errors on cTmaj are also calculated when using the 68% confidence limits for r] (grey). 
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Figure 3. Derivation of the ratio ri = ctz/ctr. In (a), the x2 statistic is used to find the best 
correspondence between the model and measured data sets for NGC 3949 (black line) and NGC 
3982 (grey line). The 68% (Icr) confidence limits on r) are determined by an increase of by 
1 (dotted lines). The disk heating of NGC 3949 (black cross) and NGC 3982 (grey cross) are 
placed on Figure 5 from Shapiro et al. (2003) in (b). 
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